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SUMMARY 

An  underwater  towed  sound  generator,  used  in  sonar  research  work,  was  found  to 
be  unstable  when  under  tow.  Wind  tunnel  tests  were  therefore  carried  out  on  a  model 
of  the  vehicle  to  assess  the  effectiveness  of  various  modifications  designed  to  improve 
its  stability. 

Fluid  dynamic  data  from  these  tests  has  been  incorporated  into  a  computer  model  of 
the  towed  vehicle  and  cable.  This  document  describes  and  presents  predictions 
obtained  from  the  computer  model  for  a  range  of  configurational  and  operational 
variables. 
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1.  INTRODUCTION 

An  underwater  towed  sound  generator,  used  in  sonar  research  work,  was  found  to  be 
unstable  when  under  tow.  A  model  of  the  vehicle  was  therefore  tested  in  a  wind 
tunnel  to  establish  the  fluid  dynamic  characteristics  of  the  vehicle  and  to  assess 
various  modifications  intended  to  increase  its  stability.  This  work  was  reported  in 
reference  1.  The  data  measured  in  the  wind  tunnel,  when  combined  with  the  other 
known  characteristics  of  the  vehicle  and  tow  cable,  is  sufficient  to  enable  predictions 
to  be  made  of  the  behaviour  of  such  vehicles  when  under  steady  tow. 

This  report  presents  a  computer  model  of  the  underwater  towed  vehicle  and  cable 
system  which  was  developed  to  use  the  wind  tunnel  measured  data.  The  model 
predicts  the  important  parameters  of  the  towed  system  for  a  range  of  conceptual 
configurations  and  operational  variables.  Results  are  presented  for  two 
configurations  which  were  shown  by  the  wind  tunnel  tests  to  be  suitably  stable. 


2.  CONFIGURATIONS  SELECTED  FOR  STUDY 

The  two  configurations  for  which  predictions  are  presented  are  shown  in  figure  1. 
They  consist  of  the  original  vehicle  body,  nose  and  strongback  fitted  with  either  of 
the  two  tail  units  shown.  (The  tail  units  are  referred  to  as  "tail  1”  and  "tail  3"  to 
match  the  nomenclature  used  in  reference  1).  The  use  of  a  cut-down  strongback  as 
indicated  in  figure  1  is  also  considered  and  is  commented  on  later. 

Details  of  the  tow  cable  configuration  are  shown  in  figure  2;  this  is  the  existing  tow 
cable  with  no  modifications.  The  arrangement  of  clip-on  discs  and  fairings  is 
designed  to  suppress  cable  strumming  (a  lateral  vibration  caused  by  vortex 
shedding),  and  they  cover  the  entire  submerged  length  of  the  cable. 


3.  THE  PREDICTION  PROGRAM 

A  listing  of  the  prediction  program  is  given  in  Appendix  1.  The  program  is  written 
in  FORTRAN  77  and  incorporates  all  input  data  and  options  within  itself.  Thus 
successive  runs  for  differing  configurations  require  program  statements  to  be  edited 
and  the  program  recompiled  between  each  run.  This  is  not  a  serious  handicap  since 
a  single  run  is  sufficient  to  provide  all  the  relevant  data  for  one  vehicle  configuration 
(see  Appendix  2  which  gives  a  sample  of  the  program  output). 

The  assumption  is  made  that  for  reasons  of  symmetry  the  underwater  vehicle  will,  if 
stable,  tow  with  negligible  sideslip  and  with  the  strongback  on  the  top.  The  program 
therefore  considers  only  forces  and  moments  in  the  vertical  plane.  The  program 
first  determines  the  incidence  at  which  the  vehicle  will  tow  by  calculating  all  relevant 
forces  (weight,  buoyancy,  fluid  dynamic  forces  and  tow  cable  tension)  and  then 
searching  for  the  incidence  at  which  all  forces  and  moments  sum  to  zero.  If  this  is 
not  achieved  within  an  incidence  range  of  -30°  to  30°  the  program  prints  "vehicle  is 
unstable".  Although  this  is  not  necessarily  true  it  does  imply  that  stability,  if  it  is 
achieved,  is  at  an  undesirably  high  angle  of  incidence  where  no  measured  data  is 
available  and  where  asymmetric  shedding  of  body  vortices  could  lead  to  dynamic 
stability  problems. 
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Once  the  incidence  of  the  towed  vehicle  is  determined  the  program  predicts  the 
tension  and  shape  of  the  tow  cable  by  considering  it  in  one  metre  lengths,  each 
length  taken  to  be  a  rigid  section  freely  pin  jointed  to  its  immediate  neighbours.  For 
each  section  of  cable  the  program  calculates  all  relevant  forces,  performs  a  force 
balance,  and  thus  determines  the  progression  of  tension  and  angle  along  the  entire 
cable. 


4.  GENERATION  OF  INPUT  DATA  FOR  THE  PROGRAM 

The  input  data  for  the  program  consist  of  all  the  parameters  of  the  tow  system  which 
are  pertinent  to  its  behaviour  under  tow.  With  the  exception  of  the  vehicle  fluid 
dynamic  characteristics  these  quantities  are  listed  in  table  1.  Where  possible  these 
quantities  have  been  directly  measured.  Otherwise  they  have  been  calculated  from 
realistic  starting  assumptions  (eg  vehicle  mass  and  buoyancy  data  are  based  on 
measured  results  for  the  existing  body  and  strongback,  plus  data  calculated  for  the 
add-on  tail  units  assuming  they  are  fabricated  from  aluminium  alloy).  The  drag 
coefficients  used  for  the  tow  cable  were  estimated  from  data  in  reference  2. 

The  fluid  dynamic  characteristics  of  the  vehicles  were  derived  from  wind  tunnel  test 
data  (see  reference  1)  and  are  input  to  the  program  as  constants  in  fourth  order 
polynomials  fitted  to  the  measured  data. 


5.  PREDICTIONS 

Ail  the  predictions  given  here  apply  only  to  a  system  in  a  steady  state.  No  attempt 
has  been  made  to  model  the  inevitable  launch  and  deployment  transients,  and 
therefore  the  program  can  give  no  information  on  whether  these  transients  will 
damp  out  and  allow  the  system  to  reach  a  steady  state.  However,  for  the  relatively 
simple  shapes  and  low  attitudes  considered  here  there  is  no  reason  to  suspect  any 
significant  dynamic  instabilities  to  occur. 

5.1  Sideslip  angle  of  vehicle 

In  the  sideslip  plane,  gravity  and  buoyancy  forces  have  no  effect  and 
therefore  only  the  fluid  dynamic  forces  need  be  considered.  Since  both 
vehicles  are  symmetric  and  statically  stable  it  follows  that  if  a  steady  state  is 
reached  (see  earlier  comment),  it  will  be  at  a  sideslip  angle  of  zero  and  with 
the  strongback  on  the  top  of  the  vehicle.  In  practice,  of  course,  sideslip 
angles  will  fluctuate  over  a  small  range,  the  size  of  the  fluctuation  depending 
on  turbulence  in  the  water,  off-axis  loads  from  the  tow  cable  and  the  size  of 
the  fluid-dynamic  stabilising  moments.  For  both  vehicles  studied  here  the 
fluid  dynamic  stability  in  the  sideslip  plane  (see  ref  1)  is  sufficient  to  ensure 
that  sideslip  angles  remain  small.  The  use  of  a  cut-down  strongback  would 
further  enhance  this  stability. 
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52  Incidence  angle  of  vehicle 

The  incidence  angle  of  the  vehicle  under  tow  will  depend  on  gravity, 
buoyancy  and  tow  cable  forces,  as  well  as  the  fluid  dynamic  forces. 
Predictions  of  this  angle  for  a  range  of  tow  speeds  and  tow  attachment  points 
are  given  in  figure  3.  Ideally,  the  vehicle  should  tow  at  a  low  angle 
throughout  the  speed  range,  and  it  suggested  that  a  limit  of  ±20°  should  be 
safely  allowable  (at  angles  of  30°  and  greater  asymmetric  growth  and 
shedding  of  body  vortices  could  introduce  some  dynamic  instability). 
Accepting  this  limitation,  figure  3  shows  that  the  vehicle  incidence  is  very 
sensitive  to  tow  attachment  point  and  that  to  achieve  acceptable  angles 
throughout  the  speed  range  requires  the  tow  attachment  to  be  within  about 
0.1m  of  the  optimum  point. 

An  optimum  tow  attachment  point  on  the  strongback  (ie  a  point  for  which 
the  incidence  stays  close  to  zero  throughout  the  speed  range)  may  be  defined 
as  the  point  for  which  the  vehicle  incidence  is  within  -2°  to  +  2°  when  fully 
submerged  and  at  rest.  Figures  3a  and  b  appear  to  indicate  that  this  can  be 
approximated  by  supporting  on  the  strongback  at  a  point  above  the  CG,  but 
figure  3c  is  included  to  show  that  this  is  not  necessarily  so.  Figure  3c  shows 
results  for  a  vehicle  fitted  with  nearly  50  kg  of  ballast  in  the  nose,  which  shifts 
the  centre  of  mass  to  about  0.13m  ahead  of  the  centre  of  buoyancy.  The 
optimum  support  point  for  this  configuration  is  a  point  nearly  0.1m  ahead  of 
the  CG. 

Reference  1  shows  that  the  use  of  a  cut-down  strongback  has  only  a  very 
small  effect  on  the  incidence  plane  fluid  dynamic  characteristics.  Figure  3 
can  therefore  also  be  used  for  corresponding  configurations  with  a  cut-down 
external  strongback,  where  the  removed  portion  is  replaced  with  an  internal 
strongback  of  roughly  equivalent  mass  and  displacement. 

5.3  Tow  Cable  Behaviour 

Factors  influencing  the  behaviour  of  the  tow  cable  are  the  loads  exerted  by 
the  towed  vehicle  plus  the  gravity,  buoyancy  and  fluid  dynamic  forces 
distributed  over  the  whole  length  of  submerged  cable.  Predictions  of  cable 
trajectory,  tension  and  total  length  are  plotted  against  depth  in  figures  4,  5 
and  6  for  varying  tow  speeds.  Note  that  the  depth  scale  reads  from  the  towed 
vehicle  (depth  =  0)  up  towards  the  surface.  The  predictions  thus  show  the 
cable  characteristics  along  the  full  length  of  submerged  cable  up  to  a 
maximum  depth  of  140m,  and  can  be  simply  truncated  at  any  intermediate 
point  for  shallower  towing  depths.  The  predictions  are  given  for  a  near 
optimum  towline  attachment  position  (which,  for  these  configurations,  is  at 
the  CG  of  the  towed  vehicle),  and  for  two  other  acceptable  positions  at  the 
highest  tow  speed  to  give  an  indication  of  the  effect  of  varying  the  tow 
attachment  position. 
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As  expected,  the  tow  cable  conditions  are  very  similar  for  both  vehicles 
considered,  and,  as  long  as  the  vehicle  remains  stable,  are  not  particularly 
sensitive  to  changes  in  tow  attachment  position.  At  low  speeds  (up  to  3  m/s) 
the  cable  tension  is  dominated  by  the  unbuoyed  weight  of  the  vehicle.  At 
higher  speeds,  and  for  deep  towing,  the  maximum  cable  tension  is 
increasingly  dependent  on  the  drag  of  the  cable  itself.  It  is  only  at  high 
speeds  (4  to  5  m/s)  and  shallow  depths  (up  to  50m)  that  the  maximum  cable 
tension  is  significantly  dependent  on  the  fluid  dynamic  characteristics  of  the 
vehicle.  It  is  therefore  reasonable  to  expect  that  for  normal  (low  speed) 
towing  conditions  the  predictions  would  be  little  changed  by  the  use  of  a  cut- 
down  strongback  on  the  towed  vehicle. 


CONCLUDING  REMARKS 

Both  of  the  vehicles  studied  here  can  be  arranged  to  tow  in  a  stable  manner 
throughout  the  desired  range  of  speeds  and  depths. 


For  optimum  behaviour  of  the  vehicle  under  tow  (ie  attitude  remaining  close 
to  zero  throughout  the  speed  range)  the  tow  cable  should  be  attached  at  a 
position  which  causes  the  vehicle  attitude  to  lie  within  the  range  of  about  -2° 
to  +2°  when  the  vehicle  is  at  rest  and  fully  submerged.  (Note  that  unless  the 
centre  of  mass  and  the  centre  of  buoyancy  of  the  vehicle  are  coincident,  then 
the  vehicle  attitude  in  air  will  not  be  the  same  as  in  water). 

If  the  tow  cable  cannot  be  attached  at  the  optimum  position  then  the  selected 
attachment  point  should  be  acceptable  if  the  vehicle  attitude  lies  within  a 
range  of  -20°  to  +  20°  when  at  rest  and  submerged. 

For  moderate  changes  in  the  vehicle  fluid  dynamic  characteristics  (eg  a  shape 
change  to  nose,  strongback  or  tail),  the  tow  cable  conditions  should  not 
change  significantly  as  long  as  the  tow  cable  is  correctly  attached  and  the 
vehicle  remains  stable  under  tow. 

The  disks  which  are  part  of  the  anti-strumming  attachments  to  the  tow  cable 
make  the  cable  a  potentially  high-drag  component  of  the  towed  system.  This 
is  not  a  problem  at  low  speeds  (up  to  3  m/s)  or  for  shallow  depths  (up  to 
50m),  but  could  become  very  significant  at  higher  speeds  and  greater  depths. 
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APPENDIX  1.  LISTING  OF  THE  PREDICTION  PROGRAM 

c  Program  name  is  TOW2.  FOR 

c 

c  This  program  predicts  the  trim  incidence  and  tow  cable  trajectory 
c  for  an  underwater  towed  low  frequency  sound  generator 
c  (with  4  cylindrical  fins  type  tail) 
c 

DIMENSION  ctens(  14),posx(  14),posy(  14),length(  14) 
c 

REAL  mass,nor0,norl,nor2,nor3,nor4,norfor,mres,length 
c 

INTEGER  dirm,dirml 
c 

common  dens,vel 

Common/hydfor/nor0,norl,nor2,nor3,nor4,ax0,axl,ax2,ax3,ax4, 

+  xcp0,xcpl,xcp2,xcp3,xcp4,refa,refl 

common/towfor/weight, bouyf, norfor,axfor>xcg,xbouy,xcp,xtow, ztow 
common/cable/tax,tnor,alph,cdcnor,cdcax,crefa,effcwt 
c 

c  set  the  fluid  constants 
c  dens  =  fluid  density(kg/ m*  *3) 

dens  =  1025.0 
c 

c  set  the  vehicle  constants 
c  refl  =  ref.  length  =  vehicle  diam(m) 
c  refa  =  ref. area = vehicle  cross  sectional  area  (m**2) 
c  mass = vehicle  mass(kg) 
c  disp = vehicle  displacement^*  *3) 
c  bouyf = bouyancy  force  on  vehicle(n) 
c  xcg= vehicle  eg  position  aft  of  nose(m) 
c  xbouy = centre  of  bouyancy  position  aft  of  nose(m) 
c  ztow  =  distance  from  vehicle  axis  to  tow  attachment  point(m) 
refl  =  0.450 

refa  =  3.1416*(refl**2)/4.0 
mass =330.3 
weight  =  9.807 *  mass 
disp =0.1584 
bouyf =9.807*dens*disp 
xcg=  1.239 
xbouy  =  1.244 
ztow =0.305 
c 

c  set  the  hydrodynamic  constants  -  these  are  the  constants  in  the 
c  following  equations: 
c  normal  force  coeff( + ve  upwards) 

c  cnor  =  nor0  +  norl’inc  +  nor2*inc**2  +  nor3*inc**3  +  nor4*inc**4 

c  axial  force  coeff( + ve  backwards) 

c  cax=  axO  +  axl’inc  +  ax2*inc**2  +  ax3*inc**3  +  ax4*inc**4 
c  centre  of  pressure  of  normal  force(calibers  aft  of  nose) 
c  xcp=  xcpO  +  xcpl’inc  +  xcp2,inc*,2  +  xcp3*inc**3  +  xcp4*inc**4 
c  WHERE:  inc=  incidence  (radians) 
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norO  =  0.030 
norl  =  12.676 
nor2  =  0.334 
nor3  =  -9.98 1 
nor4  = -2.748 
ax0  =  0.500 
ax  1= -0.078 
ax2=  1.424 
ax3  =  0.251 
ax4  =  -4.592 
xcp0= 4.232 
xcpl  =0.010 
xcp2=- 1.088 
xcp3  =  -0.139 
xcp4=- 1.153 
c 

c  set  the  cable  constants 

c  cdcax  =  drag  coeff  of  cable(per  metre)  for  flow  along  cable 
c  cdcnor=drag  coeff  of  cable(per  metre)  for  flow  normal  to  cable 
c  cdrefa=ref  area  for  cable  drag  =frontal  area  of  bare  cable/n(m**2) 
c  cmass  =  cable  mass /m( kg) 

c  cdisp= cable  displacement/m(m**3) 

c  effcwt  =  effective  cable  weight/m  after  allowing  for  bouyancy(n) 
cdcax  =  0.24 
cdcnor  =  0.16 
crefa=0.016 
cmass  =  0.763 
cdisp  =  0.000468 

effcwt  =  (cmass-cdisp*dens)*9.807 
c 

c  set  initial  values  of  variables 
c  vel  =  tow  speed(m/s) 

c  xtow  =  towline  attachment  position  aft  of  nose(m) 

vel=5.0 
xtow=xcg-0.5 
c 

c  set  up  output  headings 
write(6,2000) 

2000format(’0’,35x, ’UNDERWATER  SOUND  GENERATOR  TOWING’, 
1  ’CHARACTERISTICS’,//) 
c 

c  set  an  initial  guess  for  incidence  angle  and  incidence  step 
c  alph  =  incidence  angle(radians),  nose  up  is  +ve 
c  dalph  =  incidence  step(radians)  to  be  used  to  refine  guessed  alph 
20  alph  =  -30.0/57.296 
dalph  =  10.0/57.296 
c 

c  set  misc  indicators 

c  dirm=  1  for  +ve  pitching  moment,  =  -l  for  -ve 
c  dirml  =  the  previous  (last)  value  of  dirm 
dirm  =  1 
dirml  =  1 
c 

c  calculate  the  hydrodynamic  forces  on  the  vehicle 
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10  call  hydfor(alph,norfor,axfor,xcp) 
c 

c  calculate  towing  forces(normal  and  axial)  and  resultant  moment  on 
c  vehicle 

call  towfor(alph,tnor,tax,mres) 
c 

c  adjust  incidence  appropriately  as  determined  by  unbalanced  moment 
dirml  =  dirm 
if(mres.gt.O)then 
dirm=  1 

else 

dirm  =  -l 
end  if 

if(dirm.ne.dirml)dalph  =  -dalph/ 10.0 
c  if  close  to  stable  trim  (dalph.  lt.0.1  deg)  jump  out  of  incidence 
c  adjusting  loop 

if(abs(dalph).lt.0.001)go  to  100 
alph  =  alph  + dalph 

c  if  unstable  for  alph  between  -30  deg  and  +30  deg,  jump  out  of 
c  loop  and  write  appropriate  message 
if(abs(alph).gt.0.55)go  to  1 10 
go  to  10 
c 

c  printout  final  prediction 
c  jump  to  here  if  unstable 
1 10  write(6, 1010)vel,(xcg-xtow) 

1010  format(’  ’,’for  tow  speed  \f4.1,’  m/s  and  tow  position  \f5.2, 

1’  m  ahead  of  eg,  vehicle  is  unstable’) 
go  to  120 

c  jump  to  here  if  stable 
100  alphd  =  alph*57.296 

write(6,1000)vel,(xcg-xtow),alphd 
1000  format(’  ’,’for  tow  speed’, f4. 1/  m/s  and  tow  position  \f5.2, 

1’  m  ahead  of  eg,  vehicle  trim  angle  is  \f6.1,’  deg’) 
c  now  calculate  and  print  out  the  cable  trajectory 


call  cable(length,ctens,posx,posy) 

write(6,1020)length 

write(6,1030)ctens 

write(6,1040)posx 

write(6,1050)posy 

1020 

format(’  ’,’ 

cable  length  (m) 

\14f7.1) 

1030 

format(’  ’,’ 

cable  tension  (n) 

\14f7.0) 

1040 

format(’  ’,’ 

horizontal  position  (m)  ’,14f7.1) 

1050 

120 

c 

format(’  ’,’ 
continue 

vertical  position  (m) 

\14f7.1,/) 

c  now  set  new  tow  speed  and/or  tow  position  and  repeat  calculations, 
c  if  all  speeds  and  tow  positions  completed,  go  to  STOP 
vel=vel-1.0 
if(vel.lt.0.0)then 
write(6,*) 
write(6,*) 
write(6,*) 
vel  =  5.0 
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xtow  =  xtow  +  0.1 
if(xtow.gt.(xcg  +  0.51))go  to  200 
go  to  20 
else 

go  to  20 
end  if 
c 
c 

200  STOP 
END 


SUBROUTINE  HYDFOR(alph,norfor,axfor,xcp) 
c 

c  This  subroutine  calculated  the  hydrodynamic  forces  acting  on  the 
c  towed  body  from  curves  of  coefficients  versus  incidence 
c 

c  alph  =  incidence  of  vehicle(radians),  +ve  is  nose  up 
c  nort'or  =  normal  force  on  vehicle(n),  +  ve  is  up 
c  axfor  =  axial  force  on  vehicle(n),  +  ve  is  backward 
c  xcp  =  centre  of  pressure  of  normal  force(cal),  +ve  is  behind  nose 
c 

real  nor0,norl,nor2,nor3,nor4,norfor,norfcf 
c 

common  dens,vel 

Common/hydfor/nor0,norl,nor2,nor3,nor4,ax0,axl,ax2,ax3,ax4, 

+  xcp0,xcp  I,xcp2,xcp3,xcp4,refa,ref  1 
c 

dynp  =  0.5  *  dens  *  (vel  *  *  2) 

norfcf  =  nor0  +  nor  1  *alph  +  nor2*(alph*  *2)  +  nor3*(alph*  *3)  +  nor4*(alph*  *4) 
norfor  =  norfcf*  dynp*  re  fa 

axfcf=ax0  +  axl*alph  +  ax2*(c.lph**2)  +  ax3*(alph**3)  +  ax4*(alph**4) 
axfor  =  axfcf*dynp*refa 

xcpcal  =  xcp0+ xcpl’alph  +  xcp2*(alph*  *2)  +  xcp3*(alph*  *3)  +  xcp4*(alph*  *4) 

xcp=xcpcal*refl 

RETURN 

END 


SUBROUTINE  TOWFOR(alph,tnor,tax,mres) 
c 

c  This  subroutine  calculates  the  tow  cable  forces  on  the  towed  body 
c  and  also  the  resultant  pitching  moment  on  the  vehicle 
c 

c  alph  =  incidence  of  vehicle(radians),  +  ve  is  nose  up 
c  tnor= cable  force  normal  to  vehicle  axis(n),  +ve  is  up 

c  tax  =  cable  force  parallel  to  vehicle  axis(n),  -t-ve  is  forward 

c  mres = pitching  moment  on  vehicle(n.m),  +  ve  is  nose  up 

c 

c 


c 


common/towfor/weight,bouyf,norfor,axforpccgtxbouytxcp,xtow,ztow 
real  norfor, mres 
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tnor  =  ( weigh  t-bouyf)*cos(alph)-norfor 
tax  =  (weight-bouyf)*sin(a!ph)  +  axfor 

mres  =  (weight*xcg-bouyf*xbouy)*cos(alph)-norfor*xcp-tnor*xtow- 
ltax*ztow 
RETURN 
END 


SUBROUTINE  CABLE(length,ctens,posx,posy) 
c 

c  This  subroutine  calculates  the  tension  and  shape  of  the  first 
c  2000  metres  of  tow  cable  from  the  towed  body  by  considering  the 
c  cable  to  be  made  up  of  2000  pin-jointed  rigid  sections,  each  section 
c  one  metre  long 
c 

c  length  =  length  of  cable(m),  measured  from  vehicle  attachment  point 
c  ctens  =  cable  tension(n),  at  points  defined  by  ‘length1 
c  posx  =  horizontal  distance(m),  from  attachment  point,  +ve  is  forward 
c  posy = vertical  distance(m),  from  attachment  point,  +  ve  is  up 

c 

dimension  ctens(  14),posx(  14),posy(  14),length(  14) 
c 

common  dens,vel 

common/cable/tax, tnor, alph,cdcnor,cdcax,crefa,effcwt 
c 

real  length 

c 

c  set  up  initial  values  at  cable  attachment  point 
xold  =  0.0 
yold  =  0.0 

told = sqrt(tax*  *2  +  tnor*  *2) 
angl  =  atan(tnor/tax) 
if(tax.  lt.0.0)angl  =  angl  +  3. 14 16 
angold  =  alph  +  angl 
ystore  =  9.6 
j=0 
c 

c  now  calculate  cable  characteristics  at  each  metre  for  2000  metres 
do  10  i  =  1,2000 
xnew =xold  +  cos(angold) 
ynew = yold  +  sin(angold) 

fnor=cdcnor*crefa*0.5*dens*(vel*sin(angold))**2 
fax = cdcax  *  crefa*  0.5  *  dens  *  (vel  *  cos(angold) )  *  *  2 
wnor  =  effcwt*cos(angold) 
wax  =  effcwt  *sin(angold) 
tnewax  =  told  +  fax  +  wax 
tnewnr  +  wnor-fnor 
tnew  =  sqrt(tnewax*  *2  +  tnewnr*  *2) 
angnew  =  angold  +  atan(tnewnr/tnewax) 
c  store  variables  for  output  at  10  metre  depth  increments 
if(ynew.ge.ystore)  then 

H+l 

if(j.eq.l5)  go  to  20 
length(j)  =  i 
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10 

20 


ctens(j)  =  tnew 
posx(j)  =  xnew 
posy(j)=ynew 
ystore =ystore  + 10.0 
end  if 
xold  =xnew 
void  =ynew 
told  =  tnew 
angold  =  angnew 
continue 
RETURN 
END 
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APPENDIX  2.  SAMPLE  OUTPUT  FROM  THE  PREDICTION 

PROGRAM 


avoKawAfts  louao  amutoi  towxpa  aijumiuna 


t  at 

Caw  spaod  3.0  a/s  sad  Caw 

past tioa 

0.30 

a  aOoad 

•  t  Cf. 

vobicio  t ria 

aaf  la 

is  0 . 

.  9  dap 

csalo  Isaptb  isi 

14.3 

3*  .  0 

41 . 3 

44.3 

•  4.3 

104.3 

123.3 

144.3 

149.3 

191.3 

213.0 

231.3 

243  .3 

231  .3 

caoia  coasros  <at 

1*37. 

23*4  . 

2*47. 

3412. 

4340  . 

5131. 

5*57. 

4791  . 

7730  . 

•  411. 

9423  . 

10340  . 

11427 . 

124J7. 

sortsoatsl  position  la) 

1.7 

21 . 1 

3  4.# 

4*.# 

47 . 0 

*4.4 

102 . 9 

121.3 

142.1 

141.1 

1*1  .3 

204  .  « 

227  .1 

249  3 

vortical  pool ci jo  (at 

10 .1 

it. 7 

29.4 

J*  .  7 

50.1 

40.3 

70  .  0 

79.4 

90.3 

99.4 

109 . 9 

119.4 

120.3 

139  .  4 

far 

cow  apsad  4.3  a/s  sad  tow 

pasltiaa 

0.30 

a  aOaad 

at  Cf , 

vablcls  tria 

aaql  a 

is  1  , 

.  7  dap 

cools  UdfCtt  1  SI 

12.3 

23.3 

39.3 

54.3 

4*  .  0 

14 . 3 

103.3 

123.3 

129  .3 

13*.  3 

177.0 

194.3 

217.3 

23  7  .3 

caolo  Csasioa  :al 

i4i: . 

1434  . 

1**3. 

2111. 

24*1  . 

2*34. 

1233. 

3427. 

4074  . 

433  4  . 

49V*  . 

3449  . 

5*94  . 

4500  . 

borisoncsl  pasitioa  (a) 

4.4 

14.* 

24.4 

35.1 

0.3 

40.7 

74.3 

I*  .  4 

104.5 

123.4 

134  .* 

133  .  2 

171  .  4 

114.1 

vortical  pasltiaa  (at 

10.9 

20.1 

30  •  1 

43.2 

49.* 

43.1 

70  .1 

79  4 

90 .3 

130.1 

'.OP  .9 

119.4 

130 . 1 

129.9 

for 

cow  apaad  3.3  a/s  aad  Caw 

paiitioo 

0.30 

a  abood 

of  Cf , 

vobicio  tria 

sap  la 

is  3  . 

.  3  dsf 

caolo  1  o  o  p  t  A  at 

11.9 

22.3 

34.3 

44.3 

5  *  .  3 

71.3 

13.3 

#9 . 3 

113  .  0 

121.0 

143 .  a 

131.0 

174.3 

190.3 

caolo  taosioa  at 

1211. 

127*. 

13  45  . 

1*41. 

1344  . 

141*  . 

1*29  . 

1977  . 

2131. 

230  2  . 

247*. 

24  59  . 

2*33  . 

J0  5  3  . 

borisoQCsl  paaitloi  (at 

4.2 

*  .  2 

15.4 

22.2 

29.3 

37.* 

47.3 

37.0 

47.0 

7*  .0 

t*  .  2 

100 . 4 

112  .  9 

12  3  4 

vortical  pasitioa  (at 

l».l 

1*  .* 

30.2 

40 . 1 

49.4 

39.4 

70 . 0 

10 . 0 

19  .* 

100.0 

110.0 

119.7 

130.3 

140.3 

far 

tow  ipo«d  2.3  a/s  sod  tow 

pasitioa 

0  .  so 

a  aboad 

of  Cf . 

vobicio  trta 

aaf  la 

is  4 

■  2  daf 

caoia  looped  at 

10.3 

21 . 0 

U.O 

42.3 

53  .  0 

44 . 0 

73.3 

14.3 

97.0 

109.0 

120 . 0 

132.0 

144.3 

134.3 

caolo  Coos l«a  at 

U«4  . 

UI0. 

1214. 

123  3  . 

1294  . 

13  34  . 

13  79  . 

1  42  4  . 

1470  . 

1321. 

1549  . 

14  23  . 

1471  . 

1  73  4  . 

ioriiaatsi  pasltiaa  (at 

2.4 

5.5 

4.4 

12.4 

14  .  5 

20  .  » 

23.4 

30.3 

33.4 

.  41.3 

47.0 

53.2 

59.4 

44  .1 

vortical  poiitisa  (at 

1.7 

20  .  J 

2*  .  1 

40.1 

50.3 

40.4 

70.3 

ta .  2 

19.9 

100.4 

109.9 

123.2 

133.2 

143.4 

far  'I*  ip*«d  1.3  a/s  tad  Cow  patitloo  3. 30  •  abssd  a<  z  p ,  is  unstsbla 

for  cow  ipsad  3-3  s/a  add  Cow  position  3.30  ■  aOvad  at  cf.  vaaicla  Is  uastaoia 

T 


1 


foe  tow  it**4  4.0  •/»  *•* 

caSlo  laaftb  ISI 
CSblo  (IM1M 

laaMtctl  position  (at 
*amc«l  position  isi 


paaitisa  -0.**  ■  sfcos4  at  c* 

11.0  21.0  22.0  44.0 

ji*j.  >*»i.  mi*,  i”»- 

l.«  7.4  11. 1  11. 1 

10.4  II. «  !*•* 


for  too  spood  4.0  a/s  sa4  taw 
caaio  1  oaf cs  (at 
coals  tsasios  (a) 
Bonsoacsi  pssitloa  (a) 
tactical  yaaitlsa  l»> 


9  U.O  21.0  3  2.0  41.® 

4114.  411*.  41*1.  1*11. 

).l  4.1  11. 1  11. 1 

il.l  11.1  I*. I  **•! 


.•ftlCl*  4tt»  *S«1*  4*  “*.l  *•* 

11. g  41.1  11.1  11.*  !*»•• 

1147.  *111.  *40*.  ***»•  *1>*. 

11.1  11.1  17.4  *«•»  41.1 

41.1  11.7  71.1  71.7  11.7 

..11*11  441*  #**i*  4*  *41.1  *M 

14.9  41. »  71.1  it.*  l*«« 

1171.  1171.  1111.  »H'  Mil. 

go . a  11.1  Il.l  * « - *  «l*i 

11.1  11.1  71.1  71.1  11.7 


lit.*  111.*  147.0  1*1.0  174.0 

7171.  7*0*.  7*1*.  *1*7.  Mil. 

41.4  71.4  *1.0  Il.l  1*4.1 

too.l  S01.«  110.1  111.1  140. » 


111.0  tit. I  11*. *  110.0  1*1.0 

1410.  111*.  till-  Mil-  »»• 

1* . I  17.1  11. »  71. J  *1.7 

100.4  110.0  HO-*  HI-’  HI* 


t*t  144  4***0 
(•«  • *••* 
far  tow  afotO 
far  taw  spawn 


i . o  a/s  sad  taw 
2.0  a/s  aa«  taw 
i.O  a/s  sad  taw 
0.0  a/s  is«  taw 


passtiaa  -0.«0 
pasltiaa  -0.10 
passtiaa  -O.«0 
pasltiaa  -4.40 


a  itasO  at  cf* 
a  aboad  at  (f, 
a  itasO  sf  «f* 
a  sbaan  at  »f, 


voalsia  is  uastsals 
waasala  is  uastsal* 
waaiaia  is  nastsai* 
ostial*  is  aastatl* 


far  taw  spawn  l ■ 0  */#  san  taw 
far  taw  spawn  4.0  */*  **n  taw 
far  (is  spawn  J.o  */*  «*n  taw 
far  tow  spawn  2.0  a/s  san  taw 
far  taw  spaod  1.0  a/s  san  taw 
far  ta#  spawn  0.0  a/s  san  taw 


pasltiaa  -O.S®  •  **•*«  si  if, 

passtiaa  -0.5®  a  staan  at  «f. 
pasltiaa  -0.30  a  staan  at  ef, 


vatials  la  aastatl* 

oatisla  is  aastatl • 
watlaia  is  aaststl* 


TABLE  1.  INPUT  DATA  USED  IN  THE  PREDICTION  PROGRAM 


centre  of  mass  (eg) 
centre  of  buoyancy  (cb) 


Vehicles 


General 


Tow  cable 


with  tail 


Xcg 

1.277  m 

1.244  m 

xcb 

1.266  m 

1.239  m 

mass 

340.8  kg 

330.3  kg 

displacement 

0.162  m3 

0.158  m3 

tow  speed 

variable,  Om/s  to  5  m/s 

'tow 

variable,  0.5m  to  -0.5m 

^ow 

0.305  m 

fluid  density 

1025  kg/m3  (sea  water) 

mass 

0.763  kg/m 

displacement 

0.000468  m3/m 

drag  coeffs 

0.24  (for  axial  flow) 

reference  length  is 

0.16  (for  cross  flow) 

diam  of  bare  cable  (0.016m) 

FIGURE  1 :  GEOMETRY  OF  THE  TOWED  VEHICLES 


2 


FIGURE  3.  INCI 


•ia 


'ow  Speed 

t>)  With  tail  3 


OF  THE  TOWED  VE 


O  3 


a)  with  tail  1 


FIGURE  6.  TOW  CABLE  LENGTH 


50  100  150 

Depth  (m) 

b)  with  tail  3 


FIGURE  6.  TOW  CABLE  LENGTH 
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